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Abstract 

A concept of semiclassically concentrated solutions is formulated for the multidimen- 
sional nonlinear Schrodinger equation (NLSE) with an external field. These solutions 
are considered as multidimensional solitary waves. The center of mass of such a so- 
lution is shown to move along with the bicharacteristics of the basic symbol of the 
corresponding linear Schrodinger equation. The leading term of the asymptotic WKB- 
solution is constructed for the multidimensional NLSE. Special cases are considered 
for the standard one-dimensional NLSE and for NLSE in cylindrical coordinates. 

1 Introduction 

Soliton phenomena is an attractive field of present day research in nonlinear physics and 
mathematics. Essential ingredients in the soliton theory are the nonlinear Schrodinger 
equation (NLSE) and its variants appearing in a wide spectrum of problems. Examples 
are coupled nonlinear optics [|], ^, ^, ||] , superconductivity [||, ^] , and excitation in lattice 
systems @. 

More exactly, solitons are identified with a certain class of reflectionless solutions of the 
equations integrable via the Inverse Scattering Transform method (see, for example, (8|). 
Such equations, including NLSE, are named soliton equations. At every instant a soliton 
is localized in a restricted spatial region with its centroid moving like a particle. The 
particle-like properties of solitons are also manifested in their elastic collisions. 

Soliton equations make up a narrow class of nonlinear equations, whereas a wider set 
of nonlinear equations, being nonintegrable in the framework of the 1ST, possess soliton- 
like solutions. They are localized in some sense, propagate with small energy losses, and 
collide with a varied extent of inelasticity. These solutions are termed solitary waves 
(SWs), quasisolitons, soliton-like solutions, etc. to differentiate them from the solitons 
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in the above exact meaning. The stability of the localized form of solitons and SWs and 
their elastic collisions have led to interesting physical applications. 

It is of interest to study the influence of external fields on the soliton propagation. 
To do that it is necessary to modify the original solitary equation by introducing variable 
coefficients representing an external field potential that breaks the IST-integrability. This 
problem was considered by variational methods Q, the theory of soliton perturbations (see, 
for example, reviews |l0|, |ll]]), or by using an appropriate ansatz [||, [L^] . 

An investigation of soliton-like states is a separate problem for multidimensional mo- 
dels. The 1ST constructions are unsuitable for them except for the D4 self-dual Yang- 
Mills equations and their reductions in 1 < D < 3 (see, for example, |0| for details). 
The problem is more complicated in view of the sing ular behavior of NLSE in 2D 0]. 
Nevertheless, soliton-like solutions for NLSE on a plane can exist, as shown in [|l5| in terms 
of a suitable ansatz. 

Since the methods for constructing exact solutions for multidimensional models are 
restricted if compared to the cases of 2D models, approximate methods should be used. 
An effective approach to the problem can be developed based on the WKB-method |l(| . 

Consider the generalized nonlinear Schrodinger equation for a "matter" field ^(x, t): 

1 W 

+V(x,t) - 2r\^(x,t)\ 2 W(x,t) = 0. 

Here, x S W 1 , t £ R ; V(x, t), A(x,t) are given functions; m is a real constant, r is a 
real parameter of nonlinearity; h is Planck's constant playing the role of an asymptotic 
parameter; |^| 2 = \I'*\E f , ty* is the complex conjugate of 

In such a form equation (||) was studied in Refs. ||, 15]. Since the field of application 



of (Q) is wider, the physical meaning of the quantities entering in ([T]) may be quite different 
from the quantum-mechanical one. In particular, in the one-dimensional problem of the 
propagation of an optical pulse ||, ||| we have to assume for equation (||) x = x € M 1 with x 
being a normed temporal variable and t being a normed space coordinate along which the 
pulse propagates. The function ^(x,t) is an envelope of the pulse field. Equation ([!]) 
takes the form 

r d 1 / d 



\~ ih di + 2^{ ih dx~ + A( ' x,t) ) + y M)- 2r I^M)rj^(M) = 0- (2) 

The functions V(x, t) and A(x, t) simulate the heterogeneity of the medium. The symme- 



try of equation (g) was considered in Ref. 17], and asymptotical solutions where studied 



In the present work we have formulated the concept of semiclassically localized solutions 
for ([]]) , following the ideas of Ref. |19 ] . These solutions are the multidimensional analogues 
of the soliton, like the solutions for (l|). The particle-like properties of a solitary wave are 
described in terms of the wave centroid. The latter is shown to move along with the bi- 
characteristics of the basic symbol of the correspondent linear Schrodinger equation. We 
construct asymptotical WKB-solutions for equation (|]) and consider some examples. 
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2 Semiclassical concentraited solutions 

Soliton solutions are known to show particle properties. In classical mechanics, a particle 
is completely described by its phase orbit. Therefore, it is natural to introduce a similar 
concept for the soliton-like solutions of the nonlinear Schrodinger equation (|l]). The way 
by which to introduce the phase orbit seems to be obvious enough. It is based on the fact 
that in quantum mechanics, the first moments of a state ^f(x, t) play the role of the phase 
orbit for the quantum system. Let now *£(x, t) be a solution of NLSE (|l|). The generalized 
position operators are x{= x) and their conjugate momentum variables are p{= —ihV), 

[xk,p s ] = ift5 k ,s, k,s = l,n. 

The mean value of an operator A by the solution \& is defined as 



(i&\A\i£) 



Here, ||^|| 2 = (^|^); 

(V\A(t)\V)= I y*(x,t)A(t)y(x,t)dx, 



is a function of time for every operator A{t) and it parametrically depends on h, 

(x)=x(t,h), (p)=p{t,h). (4) 

If there exist the limits: 

lim x(t, K) = x(t), lim p(t, H) = p(t), (5) 

then x(t) and p(t) are natural to be named the phase orbit of the classical system cor- 
responding to the given solution It is obvious that both the mean values (||) and the 
limit values (||) depend on the solution \& in the general case. Hence, the choice of the 
solution <!' meets the requirement for the expressions (|5|) to be a solution of the classical 
equations of motion. By analogy with quantum mechanics (see Ref. [19|] ), we can define 
the soliton-like solutions asymptotic in h —* as follows. 

Definition 1 Let z(t) = {(x(t),p(t)), ^ t ^ T} be an arbitrary phase orbit in M? n . 
We name the solution *f/(x,t,h) of equation ([JJ as a semiclassically concentrated solution 
(SCS) of the class CE s (z(t),N) (* G CS s (z(t),N)) if: 
(i) there exist the generalized limits Q 



V l*0M,fi)| 2 xf — r \^{P,t,h)\ 2 



(ii) there exist the centered moments 



A ( a %(t,h) = (A^ ), O^k^N. 



1 By generalized limit we mean the passage to the limit standardly defined in the distribution theory 
(see, for example, Ref. feo|). 
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Here, a and j3 are multiindices, |a| + \(3\ = k, ^ k ^ N, and A^l is an operator 

with the symmetrized (Weyl) symbol A^(p, x) = (p — p(t)) a (x — x (t)) . Recall that the 
multiindex a is a vector of the form a = (a.\, . . . , a n ), where aj ^ are integer numbers. 

n n 

In addition, |a| = Yl a j an d f° r a vector £ = (Ci> • • • > d) ^ M n we suppose = E[ C; J • 
j=i ^ ^ j=i 

Note that the vectors x(t) and in Definition [l] are by no means connected with 

each other. The vector z(t) = (p(t),x(t)) is named the classical phase orbit of the system. 

Definition |l] specifies the concept of solitary waves for asymptotic solutions of NLSE (|l|). 

Theorem 1 If a solution ^ of (Tj|j is semiclassically concentrated 6 C§s(z(t), N)), 
then z(t) = (p(t), x(t)) is a solution of a classical Hamilton system with the Hamiltonian 

H cl (p,x,t) = — (p-A{x,t)Y + V{x,t). 
2m \ J 

Proof. Let A be an operator, then the Ehrenfest theorem is true for the mean value 
of A: 

In particular, for the operators p and x we have: 

j t (p) = - n m„M !<*> = £<[&,,*]>. (6) 

Using the obvious relations 

([|^(f,t)| 2 ,#]) = ([|^(x,t)| 2 ,f]) = 0, 
we have from (|6]): 

where 

1 / \ 2 
Hi = — (ihV + A(x,t)) +V(x,t). 
2m V / 

With Definition [I] of the SCS further proof coincides with a similar one for the linear 
case [19]. 

Remark. Emphasize that, as follows from the theorem, the centroid of the SCS moves 
along the bicharacteristics of the linear Schrodinger equation. 



3 Asymptotic solutions 

In Section [2] we have discussed the definition and the basic features of SCS. Here, we study 
a theoretical possibility of construction of the WKB-asymptotic semiclassically concen- 
trated solutions of equation (||) on a limited time domain < t < T with /i-independent T. 
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Taking into account the form of the one-soliton solution of NLSE (see, for example, 
let us try solution of equation ([[]) in the form 



\P = p(9, x, t, h) exp 



(8) 



Here, 9 = h cr(x, t, h) is a "fast" variable; a(x,t,h), p(9,x,t,h), and S(x,t,h) are real 
functions regular in h, that is: 

S(x, t,h) = S (x, t) + hSi (x, t) H . 

The solution (||) is assumed to be localized according to Definition (|l[). 
The derivative operators d/dt and V act on the function (||) as follows: 



- lh m = - lh Wt 



d d 
~ za >t-^Zi = -^v| e=cons t - *( Vcr ) ™; 

0=const Go 00 



where a jt = da/dt. 

Henceforth we put d/dt\ e=const = d t , V\ e=const = V,d/d9 = d e . 
Substituting (§) into (||) we find: 

exp {-ihB t + S, t - ia,td e + V - |^V 2 - ^(W)5, - ^(Va ■ V)d e 

-i A(v 2 5) - -L(W) 2 d 2 ee - i-(V5 • V) - i-(VS • Va)d e + ^-(vsf 
2m 2m m m 2m 

+i-(A-V) + i-{A-Va)d e - -(A • V5) + i^-{VA) + ^-I 2 - 2rp 2 \p = 0. 
mm m 2m 2m 

Let us gather /i-free terms in this equation and put their sum to zero. In the obtained 
equation we separate real and imaginary parts and then separate the "fast" variable 9 
from others. As a result we come to the following system of equations which determines 
the leading term of the asymptotic solution: 

a t + ^((yS-A),Va) = 0, (9) 

S, t + V + ^-(vS-X) 2 = rb(t,x), (10) 
2m \ J 

^-(Va) 2 p M + 2p 3 =rbp. (11) 
2m 

-> n 
Here, (a, b) denotes the Euclidean scalar product of the vectors: Yl a jbj', the function 

i=i 

b(t, x) appears as a "separation parameter" in separating the "fast" variable 9, and b(t, x) 
is to be determined in what follows. 

Let us look for p(9, x, t, h) in the class of functions satisfying the conditions: 

lim p(9, x, t, h) = lim p e (6, x, t, h) = 0. (12) 

8— >oo 9— >oo ' 
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Integrating equation (11) in view of (|12j) we obtain: 



' 2mr 



b — p 2 p. 



(13) 



Let us put r = x 2 > that corresponds to the existence of soliton solutions in the case 
when equation (||) is reduced to the standard NLSE (one-dimensional with V(x,t) =, 
A(x,t) = 0) ||. Then b = b 2 > and further integration of ( |l3| ) results in 

(14) 



cosh {bx-sj2m(Va)- 2 (|<r(x, t) + <n(x, i)) | 



where the function a\(x,t) appears as a "separation constant" which is to be determined 
later. 

By definition, the "fast" variable must have the structure ^a(x,t). For equation (14) 
to correspond this condition it is necessary to set 



but 



2m 



const = a. 



v (V^) 2 

Without loss of generality we can put a = 1 and then we have 

1 



P 



(Ver) 2 



2mx 2 cosh (\a{x,t) + ai(x,t)) 
Equation ( |i~0| ) takes the form 



(15) 



2m 



-\ 2 



2m 



that, together with (^), is equivalent to the single complex Hamilton- Jacobi equation 



(5 + ia\ t + V + [v(5 + ia) - ^ 



0. 



(16) 



Thus, for the leading term of the asymptotic expansion (||), 

<&(x,t,ft) = v°(x,t,h) + o(h), 

we have: 



*° = p(9, x, t, ft) exp 



-S(x,t) + S 1 (x,t) 
ft 



(17) 



where p has the form (|l~5|) and the functions a\ (x, t) , Si (x, t) are determined from successive 
approximations. The function (|l7|) can be represented in the form: 



W(x,t,h) = 2 



(Va) 



2mx 2 1 + \^ (x,t,ft)\ 2 ' 



(18) 



where 



ty (x,t, ft) = exp <j - [S(x,t) + ia(x,t) + ft(Si(x,t) + ia\{x,t))] \ . 
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It can easily to show that is an asymptotic solution of the linear Schrodinger equation: 



ihV-A(x,t)) }V (x,t) = O(h a ), 



(19) 



where a = 1. Let us write a function which satisfies equation (||) to an accuracy 
of 0(h 2 ), in the form 

<S = $°(l + (20) 

Here, is determined by expression ( |i~7| ) and the function ^ l (9, x, t) is to be determined. 
As can be seen from ([0]), it is convenient to take the variable 9 in the form 

9 = -a(x,t) + (Ti(x,t). (21) 

Let us denote Re ^(fl, x, i) = u(9,x,t), lm^ 1 (9,x,t) = v{9,x,t). Substituting (g0|) 
into ([!]) and setting to zero summands at the equal powers of h. Extract real and imag- 
inary parts in the obtained equations and find a system of equations for function (^0[), 
This system includes the Hamilton- Jacob i equation (jig) for the functions S(x,t), a(x,t) 
and also the following equations determining the functions Si (x, t) and o\ (x, t) : 

Si,t + - ( (VS - X) , VSx) - - < V<r, Va x ) + J- Aa + J- (V<7, V) log (Va) 2 = 0, (22) 
m V / m 2m 2m 



m 

a M + -<(VS-„4) ,Vai) + -(Va,V5i) 
m V / m 



1 



ffl 



1 



(VAVS-A))+ [log(Vaf] + -((VS-A),V log(Va) 



m 



(23) 



The functions u(9,x,t), v(9,x,t) are given by the expressions 
p(9, x, t)u(9, x, t) = 



2m 1 1 

^w^r i(xV)tanh ' + ^ (w ' Wi) 



(24) 



H [Act + (Vcr, V log(Va) 2 )l (sinh cosh 9-e cosh 2 

12m 



p(#, x, t)v(9, x, t) 



1 2m I Ci(x,t) 1 
x 2 (Vo-) 2 1 coshl + 4 



— (V,VS-A) 
m 



+ [d t + -([VS-A),V))log(Va) 



(e sinh 9 — cosh 0) > . 



(25) 



Here, e = sign(<r) and the function Ci(x, t) is determined by successive approximations. 
Remark. From relations (22) and ( ]23| ) it follows that ([II]) is accurate to 0(h 2 ). 



Relation (18) can be considered as a transformation connecting the asymptotic solutions 
of the nonlinear equation ([l]) with the solution of the linear equation ([l9| ) for a = 2. 

So, the leading term Vf of the asymptotic solution of equation (|l|) is completely fixed 
by expressions (|15|), (|17|), (|20|), and (|21~|), where the functions S(x,t), a(x,t), Si(x,t), and 
ui(x, t) are determined by the system of equations (16), (|22"1), and (|23|) . 
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4 Special solutions 

For a more detailed study of the above asymptotic solutions (|l7|), let us consider some 
special cases of these solutions. 



4.1 One-dimensional NLSE 

Let us put x = x G R 1 , A = d 2 /dx 2 and A = V = in (||); then equation (||) takes the 
form 



d h 2 d 2 

ih— + — TT^T + 2x 2 \^(x,t,h)\ 2 



V(x,t,h) = 0. 



dt 2m dx 2 

The leading term of the asymptotic solution (|20|) is as follows: 

i 



^(x, t) = p(x, t, h) exp 



(S(x,t) + hS l (x,t)) 



where 



p(x,t,h) 



(7+ 



I 2mx 2 cosh (j^a(x,t) + a\(x,t)) 
Equations ([u]), (^2), and (p3| ) become: 



(S + ia) tt + — (S + ia) 2 x = 0, 
zm 



'S'l.t H 5*^5 

m 



1 3 

l,s &,x<T\,x + 7; — &,xx — 0, 

m Am 



1 Q J 1 Q 1 

&l,t H &,x a 'l,x H 0><->l,a: _ ~ 

m m A 



1 , _2_ , 

m a x iti o x 



(26) 
(27) 

(28) 

(29) 
(30) 
(31) 



To construct solution ( |27| ) and (28), let us look for a special solution of equation ( |29| ) in 
the form 

S 1 = u\t + a 2 x + ipo, a = Pit + f3 2 (x - x ), (32) 
where a\, a 2 , fti, (3 2 , ¥>o> and xq are real constants. Substitution (|32]) into (^) gives 

ai = — (/3f - «!), A = -— «2/?2- 
2m m 



Substituting (32) into ( J30| ) and (pT|) , we easily obtain the following complex equation for 
the functions Si and a\\ 



(Si + «<Tl) )t H (Si + Z(7l) )a . = 0. 



m 



(33) 



Equation (|33|) is a complex wave equation the solution of which can be written as 

w(x, t) = Si(x, t) + iai(x, t) = f(x — at). (34) 
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Here, a = (a 2 +if3 2 ) I'm and /(£) are analytical functions of the complex variable ( = x — at. 
Denote (3 2 = 2rj and a 2 = 2£, then solution (p7|) takes the form 

2rj 1 



xV2m cos h 



x exp 



-f [ x - x - £t ) + lmf(x - at) 
2 



(35) 



-(2£r (r - »7 )* + + fiRe/(x - at)) 

n m 



If f{x — at) = in (^J), then (|35| ) takes the form of the exact one-soliton solution of 
the nonlinear Schrodinger equation ( p6[ ) which is reduced to the standard form || when 
h = m = 1, x 2 = 1/2. 

If /(a? — at) 0, we have an asymptotic solution similar to the one-soliton solution. 



4.2 NLSE with a separated potential 

Consider equation ( |26| ) with a potential V(x,t): 

d h 2 d 2 



ih— + 



+ 2iS\V(x,t,h)F - V(x,t) 



v(x,t,h) = o. 



(36) 



(37) 



dt 2m dx 2 

Solution (|D and flU) is determined by @), ©, and 

(5 + Mr)* + -^(s + ia) 2 , + y = o. 

zm 

For the potential V of the separated form, 

V(x,t) =vo{t) + Vl (x), 
we can easily find for the system (37), (|30|) , and (|l|) two classes of separated solutions 



determining (27) and (28). 

The first class is described by the expressions: 

S(x,t) = cit — J vo(t)dt, a(x,t) = J y/2m(ci + v\(x))dx, 

Si(x,t) = c 2 t + c 3 , ax = (3/2) log |<7 2 | + mc 2 Ja^dx + c^, 
where c\, . . . , C4 = const. The second class is presented by 

dx 
p'(x) 



S{x,t) = c^t — j vo(t)dt + c 4 , a{x,t)=c\ t — mj 



+ c 2 , 



Si(x,t) = ait + a 3 + /(x) 



(Ti = a 2 t + a 4 + a(x) 



Here, c%, . . . , C4; ai, . . . , 04 = const and the functions p(x), /(a;), o(a:) are determined in 
quadratures from the following equations: 

— ip'ix)) 2 = -iyxix) + ca) + ^Mx) + c 3 ) 2 + c 2 , 



a'(x) 



m 
p'(x) 



-£rf(x) - ±-jPix) - a 2 
pyx) 2m 



[(p'(x)] + Cim ) /' (x) = c\m a 2 p (x) — ma\[p \x)\ — c\mp ix)p" {x) 
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The above solutions show the potentialities of separation of variables as applied to the 



solution of the general system of equations (16), (22), and 



4.3 Cylindrical coordinates 

The self-focusing effect of a beam of high power optical radiation propagating along the 



z-axis is described by a nonlinear Schrodinger equation of the form [22]: 

V(x,y,t,K) = 0. 



d h 2 

ih di + 2m" A± + 2 ^( x ^y^M 2 



(38) 



Here, t is the time coordinate in the coordinate system moving with a group velocity 
along the direction of the radiation propagation. Aj_ = d 2 /dx 2 + d 2 /dy 2 is the Laplace 
operator in a plane orthogonal to the radiation direction. A function V(t, x, y) simulates 
some nonstationarity and heterogeneity of the medium where the pulse is propagating. In 
the stationary statement of the problem, the variable z takes the place of t. 

Let us find asymptotic solutions of the form ( p0|) for equation ( |38l) in the cylindrical 
coordinates 



x = r cos cp, y = r sin ip. 
The leading term of the asymptotic solution ( |20| ) takes the form 



®(r, t, h) = p(r, (p, t, h) exp <J j[S{r, (p, t) + HSi(r, tp, t)] } , 



where 



p(r,<p,t,h) 



(V±a) 



2mx 2 cosh (^cr(r, ip, t) + a\(r, tp, t)) 



(39) 



(40) 



(41) 



V± is the gradient operator with respect to the variables r and ip 

The special solution, illustrating a transverse heterogeneity of the optical pulse, can be 
written as follows: 



■$(r,<p,t,K) 



■. cosh 1 



W 2m 



x exp 



UK 



+ 



ci \ , 1 . ai 

— +a 2 Jr + cibit + - log r + — + o 3 



C2 



[ l\2mh m 
where c\, 02, C3, ai, 03, b\ = const; c\ 7^ 0. 



Jt — mb\r + — + c 3 



5 Conclusion 

An outcome of the work are the asymptotic solutions of the solitary wave type for NLSE © 
obtained for a finite time interval < t < T. The question of the validity of long-time 
asymptotics requires a special consideration, since, even for the linear case, this problem 
has received rather much attention in the literature P3l. 
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Solitary wave solutions are considered here similar to quantum wavepackets. This per- 
mits one to investigate the particle-like properties of the SWs in terms of the Definition p] 
based on the Ehrenfest theorem. 

The equations of motion for the centroid of an SW are found to be independent of 
the nonlinearity factor r in (|l|), while the solution itself depends essentially on r. In such 
an approach, the centroid of the solitary wave solution moves like a classical particle in 
the external field described by the potentials (V,A) in (jl]). This completely corresponds 



to the well-known soliton properties for the one-dimensional case [12[. Such a situation 
takes place only when the nonlinearity factor is constant. If r = r(x, t), the field t, h) 
results in the appearance of additional classical variables, and the correspondent equations 
of motion become much more complicated. 

The problems discussed are beyond the scope of this work and requires a special inves- 
tigation. 
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